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ROSAT-HRI investigation of the NGC 507 X-ray halo 
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ABSTRACT 

We present an X-ray investigation of the elliptical galaxy NGC 507 in the 
Pisces cluster. We make use of archival Rosat HRI and Chandra data, and 
of previously published PSPC data, to connect the large-scale structure of the 
halo to the core morphology. Our analysis, based on a bidimensional double Beta 
model of the halo surface brightness, shows that the halo core (r < 2—3r e ) and the 
external halo (r > 3r e ) are characterised by different dynamical properties and 
suggests a different origin of the two components. The halo core has a complex 
morphology with a main X-ray emission peak, coincident with the center of the 
optical galaxy, and several secondary peaks. The spatial and spectral analysis 
of the central peak shows that this feature is produced by denser hot gas in the 
galaxy core. While both homogeneous and inhomogeneous cooling flow models 
predict a deposited mass exceeding the observed amount, our data support the 
scenario where the gas is kinetically heated by stellar mass losses. Comparison 
with previously published studies suggest that the core of X-ray extended galaxies 
is associated to the stellar distribution and has similar properties to the X-ray 
halo of compact galaxies. The secondary peaks are due instead to interactions 
between the radio-emitting plasma and the surrounding ISM, producing density 
fluctuations in the hot gas. We found that the energy input by the central radio 
source in the ISM is large enough to prevent gas cooling and may explain the 
failure of the cooling flow models. 

The total mass profile derived from the bidimensional model shows that a 
significant amount of dark matter is present at large radii. The dark halo extends 
on cluster scales and is likely associated with the whole cluster rather than with 
NGC 507. This structure is typical of many X-ray bright Early-Type galaxies 
and may explain the spatial and spectral differences with X-ray compact galaxies 


1 Universita di Palermo - DSFA - Sez.di Astronomia, P.zza del Parlamento 1, 90134 Palermo; pao- 
lillo@astropa.unipa.it, peres@astropa.unipa.it 

2 Osservatorio Astronomico di Capodimonte, Via Moiariello 16, 80131, Napoli - ITALY 

3 Harvard-Smithsonian Center for Astrophysics, High Energy Division, 60 Garden St., Cambridge, MA 
02138; pepi@head-cfa.harvard.edu, khn@head-cfa.harvard.edu 



2 


largely debated in literature. The large-scale surface brightness distribution is 
irregular and more extended in the North-East direction. The displacement of 
the cluster halo from the optical galaxy and the filamentary structures observed 
in the halo core further suggest that the galaxy may be slowly moving within the 
group potential. 

Finally we found that ~ 20% of the sources detected by Kim & Fabbiano 
(1995) in the NGC 507 halo are due to point sources while the nature of the 
remaining population is not clear. If associated with NGC 507, they can be 
either accreting binaries hosting a massive black hole or density clumps of the 
X-ray halo. 

Subject headings: X-rays: individual (NGC 507)—galaxies: halos—galaxies: 
clusters: individual (Pisces)—galaxies: jets—radio continuum: galaxies 


1. INTRODUCTION 

NGC 507 is the brightest member of a group of galaxies belonging to a poor system: 
the Pisces cluster. The latter in turn is part of the Perseus-Pisces supercluster (e.g. Wegner, 
Haynes & Giovanelli 1993; Sakai, Giovanelli & Wegner 1994), located at 66 Mpc 4 from 
our galaxy (Huchra, Vogelcy & Gellcr 1999). NGC 507 is classified as E/S0 galaxy (de 
Vaucouleurs et al. 1991). Gonzalez-Serrano & Carballo (2000) showed that, while the surface 
brightness profile follows a de Vaucouleur law out to ~ 15” from the galaxy center, at larger 
radii there is significant excess emission over the r 1 / 4 profile. Barton, de Carvalho & Geller 
(1998) studied the galaxy distribution in the surrounding of NGC 507 finding that it is part 
of a group including ~ 70 members. 

The first studies in the X-ray band based on Einstein data (Kim, Fabbiano & Trinchieri 
1992a,b) measured a luminosity comparable to those of poor clusters (Lx — 10 43 ergs s _1 ) 
with a temperature kT > 1.5 keV at the 90% level, and suggested the presence of a cooler 
core. These findings where confirmed by the Position Sensitive Proportional Counter (PSPC) 
on the ROSAT satellite which revealed extended X-ray emission out to a radius of at least 
16 arcmin. Kim & Fabbiano (1995) were able to resolve the cooler region within the central 
150”, thus supporting the presence of a massive cooling flow, and suggested that cooling 
clumps may be distributed at large radii. They estimated a hot gas temperature ~ 1.1 
keV with no sign of excess absorption over the galactic value. The data constrained the 


4 Assuming Hq = 75 km s 1 Mpc 1 
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metallicity to be not higher than the solar one. Matsumoto et al. (1997) exploited the higher 
spectral resolution and larger energy range of the ASCA satellite to separate two distinct 
components in the NGC 507 X-ray spectrum: a soft thermal emission with kT ~ 1 keV and 
a harder component with kT ~ 10 keV probably due to unresolved discrete sources. The 
ASCA data seem to favour a subsolar metallicity of ~ 0.3 Z & . The different metal abundance 
obtained by ROSAT and ASCA is explained in part by the new analysis performed by Buote 
(2000a) on the ROSAT PSPC data, which found a steep gradient in the metallicity profile 
of the hot halo. 

NGC 507 is a radio galaxy with a steep radio spectrum and weak core (Colla 1975; 
de Ruiter et al. 1986; Parma et al. 1986). It may be a source with particularly weak jets, 
or possibly a remnant of a radio galaxy whose nuclear engine is almost inactive and whose 
luminosity has decreased due to synchrotron or adiabatic losses (Fanti et al. 1987). Canosa 
et al. (1999) used the ROSAT HRI data to study the correlation between the radio and 
X-ray emission of a possible nuclear source, obtaining an upper limit for the nuclear X-ray 
flux which exceeds the expectation based on the radio data by more than two orders of 
magnitude. 

In this work we use data collected by the ROSAT HRI to study the structure of the 
NGC 507 hot halo. Our aim is to relate the large scale properties of the gaseous halo to the 
small scale structure of the nuclear region. To this end the HRI is an ideal instrument since 
it combines a large field of view with a good angular resolution. We complete our analysis 
using the recent Chandra ACIS-S observation to confirm the ROSAT results and to perform 
a spectral analysis of the halo core. Our results are compared to those obtained by other 
authors on several X-ray bright galaxies in the attempt to explain the different properties of 
X-ray extended and compact galaxies. 

Throughout this work we adopt H 0 = 75 km s _1 Mpc -1 . 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. The Data 

The NGC 507 field, including NGC 499, was observed on January 1995 with the ROSAT 
HRI (David et al. 1996) for a total time of ~ 28 ks (Table 1). The data, consisting of 2 
Observing Intervals (OBI), were reduced with the SASS7_5 versions of the ROSAT standard 
analysis software (SASS). For our data analysis, we used the IRAF/XRAY and CIAO pack¬ 
ages developed at the Smithsonian Astrophysical Observatory and at the Chandra X-ray 
Center (CXC), and other specific software as mentioned in the text. 
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To correct for the errors in the aspect time behaviour (Harris 1999), which affects 
the HR1 data processed with the SASS versions prior to SASS7_B of March 1999 (as it 
happens for our data), we run the correction routine ASPTIME (F. Primini 2000, private 
communication). 

We further tried to correct for problems related to the spacecraft wobble (Harris et 
al. 1998); however the pointlike sources detected in our held (§ 2.7) are not bright enough 
to allow an accurate determination of the X-ray centroid (crucial for the correction), when 
divided in time bins. Since the global PRF is ~ 7” FWHM we decided to use the uncorrected 
data. The HRI held of view (FOV) is shown in Figure 1. 

We have run the software developed by Snowden et al. (1994, hereafter SMB), to 
produce an exposure map and a background map for our observation. We then corrected 
for exposure time and quantum efficiency variations across the detector by dividing the raw 
data by this exposure map. 


2.2. Surface Brightness Distribution 

To study the surface brightness distribution of NGC 507 the exposure corrected data 
were rebinned in 5 x 5 arcsec pixels and then convolved with a gaussian hlter of a = 15” 
(Figure 2). The resulting image shows the presence of an extended X-ray halo surrounding 
the galaxy. The central part of the halo (the one enclosed in the dashed box) is relatively 
symmetric but elongated in the NE-SW direction. Two filamentary structures seem to depart 
from the galactic center toward NW and SE respectively. The morphology of the outer halo 
(r > 200”) is very irregular due to the low S/N ratio of the HRI image. However there 
appears to be more emission on the North and East sides with respect to the South and 
West ones. This is confirmed by the radial profiles discussed in § 2.5. In contrast with NGC 
507, the NGC 499 halo is less extended and much more symmetric. The positions of the 
sources detected by the wavelets algorithm (§ 2.7) and discussed in § 2.7, are marked with 
crosses. 

The detection of three distinct emission peaks (white crosses) in the core of the X-ray 
halo (dashed box in Figure 2) suggests a complex structure of the nuclear region. Examining 
this area in more detail (left panel of Figure 3), a number of interesting features can be 
detected: 1) a strong emission peak can be observed in the center of the X-ray halo; this 
peak is slightly elongated due to the presence of two ‘tails’ extending in the NW and the 
SE direction; 2) a secondary peak is visible ~ 50” West of the primary peak (source No.9); 
3) a region of low X-ray brightness is interposed between the primary and secondary peak; 
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4) a third peak is visible ~ 40” East of the galaxy center (source No.5); 5) additional 
diffuse emission is present on the Southern side of the nucleus, making the surface brightness 
gradient flatter than on the Northern side. 

The main peak is almost coincident with the center of the optical galaxy (Figure 3, right 
panel). The second and third peak, instead, are not related to any optical feature, although 
being within the stellar body of the galaxy. There is also some faint emission coincident with 
the position of the elliptical galaxy NGC 508, located 1.5 arcmin North of NGC 507: this 
detection has a low statistical significance (~ 2 a) but it is confirmed by the Chandra data 
analysis performed in § 2.8. 

To reveal the presence of extended low S/N features we adaptively smoothed the HRI 
image with the csmooth algorithm, contained in the CIAO CXC package, which convolves 
the data with a gaussian of variable width (depending on the local signal to noise range 
of the image) so to enhance both small and large scale structures. In the resulting image 
(Figure 4), in addition to the three emission peaks mentioned before, we can distinguish an 
elongated tail extending on the North and bending Westwards out to 200” (63 kpc 5 ), and 
an X-ray tongue protruding by ~ 160” (51 kpc) on the South-East side. These two features 
correspond to the filamentary structures visible also within the dashed box in Figure 2. 


2.3. X-ray/Radio Comparison 

NGC 507 is known to be a weak radio source. In Figure 5 we superimpose the radio 
contours at 20 cm (Parma et al. 1986; de Ruiter et al. 1986) on the smoothed X-ray image. 

The nuclear radio source, visible in the high resolution map (right panel of Figure 5) 
is coincident with the position of the central X-ray peak. However it seems to be slightly 
displaced (~ 4”) toward the South-East edge of the peak. The radio jets/lobes appear 
aligned with the direction of the secondary X-ray peaks. The western jet appears to be 
collimated up to 10-15 arcsec (~ 4 kpc) and then expands in a large lobe. The latter is 
coincident with the region of low surface brightness described in § 2.2. The higher resolution 
radio contours (Figure 5, right panel) reveal a very good agreement between the morphology 
of the Western lobe and the shape of the ‘cavity’. The secondary X-ray peak lies at the end 
of the radio lobe. 

The correlation between the X-ray and radio surface brightness for the Eastern lobe, 
which is fainter than the Western one (27 mJy compared to 40 mJy at 1.4 GHz; de Ruiter et 


5 At the the assumed distance of NGC 507, l’=19.1 kpc 



al. 1986) is only visible when we look at the high resolution map: the radio emission seems 
displaced South of the X-ray peak and falls in a region of low X-ray emission. 


2.4. Bidimensional Halo Model 


To study the complex structure of the inner halo of NGC 507 we built a bidimensional 
model using the Sherpa fitting software contained in the CIAO CXC package. The model 
was composed by adding together two bidimensional Beta components of the form: 



(1 + (r/r 0 ) 2 ) 3 ^- 0 - 5 

(1) 

where 

fl — -1- ?/ 2 

r ( x y'j — new V W 1 tfnew 

( 2 ) 

and 

Xnew = (x- x 0 ) cos(9) + (y - y 0 ) sin(6>) 

( 3 ) 


Unew = (y - Vo) COS (9) - (x - x 0 ) sin (6). 

( 4 ) 


where tjo are the coordinates of the X-ray centroid, /3 is the slope of the surface brightness 
profile, ro is the core radius, e is the ellipticity and 9 is the position angle (P.A.). The two 
components represent respectively the central X-ray peak and the extended X-ray halo. The 
2D model was convolved with the HRI PRF (David et al. 1996) and then fitted on the 
exposure corrected HRI image. To model simultaneously small and large scale components, 
we fitted the extended component on a 30 x 30 arcsec pixel image and the central component 
on a higher resolution 5x5 arcsec pixel image. 

Table 2 shows the best-fit parameters of the bidimensional model. For the large scale, 
fit the ellipticity and the P.A. of the two components were fixed to 0 because of the poor 
S/N ratio of the HRI image at large radii. For the central component fit leaving these two 
parameters free to vary yields e = 0.56 ± 0.02, 9 = 137°, and results in minor changes on the 
slope (/3 = 0.63 ± 0.01) and the core radius (r 0 = 0.32 ± 0.02). The central component is 
thus elongated in the NW-SE direction, as expected from the presence of the nuclear ‘tails’ 
(§ 2.2). The centroid of the extended component is displaced by ~ 22” (7 kpc) South-West 
of the central emission peak (Figure 6, left panel); this is in agreement with the position of 
source No. 10 detected by the wavelets algorithm (left panel of Figure 3) which corresponds to 
the centroid of the large-scale X-ray distribution. The best-values r o and (3 of the extended 
component are consistent within 1 o with those found by Kim & Fabbiano (1995); however 
our core radius is slightly larger due to the fact that we have separated the contribution of 
the central peak. 
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In the right panel of Figure 6 we show the residuals of the bidimensional model. Most 
of the structure already discussed in § 2.2 is clearly visible, in particular we can distinguish 
the secondary peaks and the X-ray excess, located on the NW side of the central peak, 
corresponding to the nuclear ‘tail’. On the Northern side of the halo there is an elongated 
X-ray feature corresponding to the Northern tail already seen in the adaptively smoothed 
image of Figure 4. 


2.5. Radial Brightness Profiles 

To study quantitatively the hot halo of NGC 507 we derived radial profiles of the X-ray 
surface brightness. We binned the X-ray counts in 2 arcsec annuli, centered on the main 
X-ray peak, up to 20 arcsec and in 10 arcsec annuli at larger radii. This allows to both 
exploit the HRI resolution in the inner regions, where the X-ray emission is stronger, and to 
have a good S/N ratio in the outer regions. The brightness profile is shown in Figure 7 as a 
thin continuous line. 

To derive the instrumental background level we used the background map produced by 
the SMB software (dot-dashed line in Figure 7). Since the latter may be a poor representation 
of the actual background when an extended source is present in the field (see for instance 
Paolillo et ah 2002), we compared the HRI data with the PSPC profile obtained by Kim 
& Fabbiano (1995). The HRI profile was rebinned in 30 arcsec annuli to match the PSPC 
resolution and the SMB background was then rescaled to obtain the best agreement between 
the two profiles in the 100-500 arcsec range. The final adopted background (dashed line) 
is in close agreement with the flattening level of the HRI profile (i.e. the brightness level 
measured at radii > 500 arcsec), thus indicating that the measured background is dominated 
by the instrumental contribution expected when no source is present. Using this background 
determination we derived the background-subtracted profile shown as filled circles. 

The inner 10 arcsec are dominated by the emission due to the central X-ray peak shown 
in left panel of Figure 3. Even though a comparison with the HRI PRF suggests that the 
peak is extended, taking into account the PRF uncertainties due to residual aspect problems 
(shaded region in Figure 7; David et al. 1996) and the asymmetries discussed below, we 
find that the central peak is marginally consistent with the presence of a pointlike nuclear 
source. However the analysis of the Chandra data (§ 2.8) performed by Forman et al. (2001a) 
unambiguously show that the central component is extended. 

Past 20 arcsec (~ 6 kpc) the X-ray profile flattens due to the presence of the extended 
component, and then decreases with a power-law profile for r > 60”. 



Figure 8 shows the model profile (§ 2.4) - convolved with the HRI PRF - superimposed 
on the HRI data. The contribution of the two components are shown by the dashed and 
the dot-dashed line for the central and extended component respectively. The PSPC profile 
measured by Kim & Fabbiano (1995) is in very good agreement with the HRI data and shows 
that the halo brightness decreases as a r -L88 power-law even outside the range covered by 
the HRI points. For comparison we also show the optical V band surface brightness profile 
measured by Gonzalez-Serrano & Carballo (2000). We can see that in two regions, for 
10” < r < 20” and 40” < r < 100”, the HRI counts are in slight excess with respect to those 
predicted by the bidimensional model. These excesses are due respectively to the presence 
of the nuclear tails discussed § 2.2 and to the X-ray emission at the edge of the SW radio 
lobe. 


We have shown in previous sections that the cluster component is perturbed by the 
radio emitting plasma and is displaced from the central peak. To evaluate the effect of these 
asymmetries on the extracted profile we compare, in Figure 9, the average (0° — 360°) profile 
with the one extracted in the NE (0° — 90°) quadrant, i.e. in the opposite direction of the 
cluster component centroid and where the influence of radio lobes is smaller. The comparison 
shows that the cluster component offset is responsible for the flattening of the X-ray profile 
in the range 10” < r < 30”. For r > 80” the X-ray profile is almost unchanged since the 
displacement is not relevant on large scales. We notice that this displacement affects the 
total profile but not those of the individual components. The consequences of this result will 
be discussed further in § 2.6. 

In § 2.2 we noticed that at large radii (r > 200”) there seems to be more X-ray emission 
on the North and East side of NGC 507 than on the South and West ones. To verify 
this fact we derived X-ray surface brightness profiles in 2 regions: a North-East sector 
(20° <P.A.< 110°) and a South-West one (110° <P.A.< 200°), and extracted counts in 30” 
annuli centered on the X-ray centroid. The resulting profiles are shown in Figure 10. Past 
200” the count rate in the North-East pie is systematically higher than the one in the SW 
pic, out to ~ 450”. This excess is significant at the 3a level. For r > 500” instead, the 
counts in the two pies are consistent within the errors. 

To check if the observed excess in the NE sector is preferentially located within NGC 
507 and NGC 499 we further divided the NE sector in two 90° pies. We found that the 
behaviour of the radial profiles is essentially the same in the two subsections in which the 
NE pie was divided, thus indicating that the excess emission is distributed through all the 
NE region. 
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2.6. Density, Cooling Time and Mass Profiles 

To derive the properties of the X-ray emitting gas we followed the method described 
by Kriss et ah (1983), deprojecting the measured emission in concentric shells. We use 
the model profile shown in Figure 8 , since the bidimensional model provides a satisfactory 
representation of the mean X-ray emission, as shown in § 2.5. This approach does not 
account for the fluctuations in the X-ray surface brightness due to the interaction with 
the radio-emitting plasma, which are perturbing the X-ray halo within 1 arcmin and for 
the displacement of the central and extended components. However, since we are mainly 
interested in obtaining the average gas properties, we will use this approximation and point- 
out, in the following discussion, the effects due to the lack of smoothness and spherical 
symmetry. 

We assume an isothermal temperature profile with kT =1.1 keV and the cooling func¬ 
tion given by Sarazin & White (1987). Moreover, since the temperature profile measured 
by Kim & Fabbiano (1995) shows the presence of a central temperature drop, we also used 
a lower temperature for the central component, obtaining an almost identical result. The 
deprojected electron density is shown in Figure 11. 

Using this density profile we calculated the cooling time r c oc nkT / n e nn A as a function 
of radius shown in Figure 12, where n e is the electron density, Uh = n e /1.21 (Sarazin 1988) 
and n = n e + rin are the hydrogen and total density and A is the cooling function. Our 
result is in good agreement with the Kim & Fabbiano (1995) estimate at large galactocentric 
distances (r > 60”) with minor differences due to their use of a non-isothermal temperature 
profile. At smaller radii our HRI data are able to sample the characteristics of the central 
component, whose cooling time is more than one order of magnitude smaller than that of 
the extended component. Both components, however, have cooling times smaller than the 
Hubble time (10 10 yr) within the cooling radius ~ 250” (80 kpc). Figure 13 (left panel) 
shows the integrated gaseous mass as a function of radius. The gaseous mass, calculated 
from the gas density profile assuming cosmic abundances (Sarazin 1988), yields ~ 10 12 M 0 
within 1000” (320 kpc). The contribution of the central component, which dominates the 
inner 15” (5 kpc), is estimated to be ~ 10 10 M 0 at r = 1000”, representing ~ 1% of the 
total gaseous mass. In Figure 13 (left panel) we also plot the amount of mass expected to be 
injected by stellar mass losses in 10 10 yrs using a mass-loss rate of 0.0078 x 10 -9 M 0 /L 0)B 
(Athey et al. 2002) and the optical profile shown in Figure 8. 

The dependence of these results on the assumption of spherical symmetry was evaluated 
comparing the average gas properties to those derived using only the NE quadrant profile 
shown in Figure 9. The main differences (see crosses in Figures 11, 12 and 13, left panel) are 
seen in the 20-100 arcsec region where, using the NE profile, the ‘shoulder’ corresponding 
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to the core radius of the extended component is smoother, thus decreasing the gas density 
and increasing the cooling time by ~ 30 — 40%. At larger radii the gas properties are almost 
unchanged since the displacement is not relevant on large scales. 


The total mass profile, shown in Figure 13 (right panel), was calculated using the equa¬ 
tion (Fabricant, Lecar & Gorenstein 1980): 


M (< r) 


rkT(r) /dlog p g dlogT\ 
Hm p G \ d log r dlog r J 


( 5 ) 


where M(< r ) is the mass contained within r, T is the gas temperature, and d log p g j d log r 
and dlogT/dlogr are respectively the logarithmic gradients of the gas density and temper¬ 
ature. 


We assumed three different temperature profiles: an isothermal profile at 1.1 keV, a 
linear fit kT = —2.3 x 10~ 4 r + 1.21 keV (r is expressed in arcsec) to account for the 
temperature decline after ~ 150”, and a power law kT = 0.04 r°' 5 + 0.7 keV to fit the central 
temperature drop. Figure 13 (right panel) shows that the total mass profiles are weakly 
dependent on the assumed temperature profile, except in the inner 100” (32 kpc) where the 
central mass obtained accounting for the low central temperature is ~ 50% smaller. The 
shaded region represents the total stellar mass calculated from the surface brightness profile 
measured by Gonzalez-Serrano & Carballo (2000) and assuming a M/L ratio ranging from 
6 to 8 M 0 /L 0 . 

For r > 50” our mass profile is in good agreement with Kim & Fabbiano (1995) and 
suggests the presence of dark matter, since it exceeds the stellar mass estimate extrapolated 
from the optical profile measured by Gonzalez-Serrano & Carballo (2000). For r< 50” the X- 
ray cumulative mass estimate yields an unphysical result: the total mass estimate is smaller 
than the stellar mass and the profile has a negative gradient for r > 10” reaching a minimum 
at ~ 20”. This behaviour is due to the flattening of the density profile which produces a 
drop in the density derivative contained in equation (5). In fact, using the profile extracted 
in the NE quadrant (dashed line in Figure 13, right panel), the mass prohle is almost in 
agreement with the stellar prohle. Thus the lack of azimuthal symmetry around the central 
X-ray peak, produced by the offset between the central and extended components and by the 
perturbed gas distribution, suggest that equation (5) is no longer valid on scales smaller than 
50” because the gas is not in hydrostatic equilibrium. The difficulties involved in deriving 
accurate mass profiles in the halo core from X-ray measurements have been already noticed 
in previous works (see for instance Brighenti & Mathews 1997) and can be explained by the 
presence of non-gravitational effects such as those discussed in § 3.3 6 . 


6 Within a few arcsec from the main emission peak we must also consider the influence of the PRF. Even 



11 


2.7. Discrete Sources 

To study the population of discrete sources in the HRI held we used the wavelets algo¬ 
rithm developed by Damiani et al. (1997a,b). Using a S/N threshold of 4.2 (corresponding 
to a contamination of 1 spurious source per held; F.Damiani, private communication) the 
algorithm detected 11 sources shown in Figure 14. The properties of these sources, as mea¬ 
sured by the wavelets algorithm, are reported in Table 3. Sources No. 10 and 11 are missing 
from the table since they are NGC 507 and NGC 499. We also notice that sources No.5 and 
9 correspond to the secondary central peaks discussed in § 2.2. 

To derive the X-ray huxes we used the PIMMS software to calculate a conversion factor 
between the HRI count rate and the incident hux. We assumed a power law spectral model 
with photon index 1.96 ± 0.11 (/hr/ 0 '™ in Table 3) following the results of Hasinger et al. 
(1993) for faint sources, and a thermal spectrum with kT = 0.52 and abundance Z ~ 0.2 Z Q 
(Ihri i n Table 3) as found by Kim & Fabbiano (1995). Assuming a galactic absorption 
of 5.3 x 10 20 cm 2 we obtain a conversion factor of 1 cps=6.925 x 10 11 ergs s 1 and 1 
cps=4.076 x 10 -11 ergs s^ 1 respectively in the 0.1-2.4 keV energy range. 

In Table 4 we cross-correlate the sources detected in the HRI held with those found 
by Kim & Fabbiano (1995). The huxes that they derive from PSPC data are reported in 
column fp§pc■ We searched for long-term variability of the matching sources comparing these 
huxes with those measured by the wavelets algorithm (Table 3). We mark as variable those 
sources whose hux varied by more than 3 a between the PSPC and HRI observations, i.e. 
between 1993 and 1995. We found that only source No.3 show signihcant variability in this 
interval. To check if the different procedures, extraction radii and background subtraction, 
significantly affect the PSPC count rate determination, we extracted from the GALPIPE 
database (Maekie et al. 1996, and references therein), which has produced a list of all 
sources present in the ROSAT PSPC archive 7 , the count rates measured by the wavelets 
algorithm on the second of the PSPC observations (RP600254a01) used by Kim & Fabbiano 
(1995). The results are in good agreement with the Kim & Fabbiano (1995) determinations, 
except for source No.4. In this case however the proximity to the NGC 499 halo and the large 
PSPC PRF may compromise the wavelets accuracy. We must notice that sources No.5 and 
9, which are likely to be produced by the interaction of the radio lobes with the surrounding 
ISM (see § 3.3), have no PSPC counterpart since they are undetected due to the lower spatial 
resolution of this instrument. 


if the bidimensional model accounts for instrumental effects, an accurate deconvolution is very difficult. 

'An online version of the GALPIPE database and documentation can be found at the D.I.A.N.A. home- 
page at the Palermo Observatory: http://dbms.astropa.unipa.it/ 
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The optical DSS image of the NGC 507/499 region was visually inspected to find optical 
counterparts of the X-ray sources within a 2a radius (Table 3) from the X-ray centroid. The 
results are summarised in column “Opt.id.” of Table 4. We also cross-correlated our source 
list with the catalogs contained in the NED database, finding that source No.7 is likely to 
be the X-ray counterpart of the radio source NVSS J012345+332554. 

Many of the sources detected in the PSPC data by Kim & Fabbiano (1995) are not 
detected in the HRI field. We derived upper limits on the HRI flux to find whether this 
result is clue to the lower HRI S/N ratio or by the temporal variability of the sources. For 
all PSPC sources not detected in the HRI field, we extracted count rates using the radii 
reported in Table 2 of Kim & Fabbiano (1995) rescaled by a factor 1/3 to take into account 
the difference between the HRI and PSPC PRF. The background was measured locally 
in an annulus of radius twice the one used for the source. The resulting fluxes (obtained 
with the same conversion factors discussed above) and 3a upper limits, are shown in Table 
5. The comparison with the PSPC fluxes shows that all sources except one are below the 
HRI sensitivity. Source No.6 is the only one which must be significantly variable since the 
expected flux is higher than the 3a upper limit indicated by the HRI data. 

The implications of these results on the nature of the discrete sources found in the NGC 
507/499 field will be discussed in § 3.4. 


2.8. The Chandra data 

Chandra observed NGC 507 on 2000 October 11 (ObsID 317) with the ACIS-S detector 
for ~ 28 ks. The Chandra field of view is smaller than the one of the Rosat HRI and is not 
azimuthally symmetric clue to the position of the ACIS ccd’s on the focal plane. Thus the 
Chandra data are not as good as the Rosat one’s to relate the core region to the large scale 
structure of the X-ray halo; however the higher angular resolution allows a more detailed 
study of the NGC 507 core. In this work we use the Chandra data to confirm our HRI results 
and to study the properties of the central peak; a complete analysis of the Chandra data is 
already in progress (Forman et al. 2001a,b). 

The adaptively smoothed and exposure corrected image of the ACIS-S3 chip is shown 
in Figure 15. The Chandra data confirm the complex morphology shown by the Rosat HRI: 
the galaxy has an extended halo with two X-ray peaks and a small gas tongue extending 
in the South-East direction, as discussed in § 2.2. The presence of the Northern ‘tail’ seen 
in Figure 4 is less evident; however the surface brightness on the Northern side extends in 
the North-West direction and its gradient is shallower than toward the South. This is also 
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in agreement with the presence of sharp Southern and Eastern edges in the gas distribution 
reported by Forman et al. (2001b), even though in the adaptively smoothed image they are 
not as evident as in the raw data. As discussed in § 2.5, the surface brightness profile of the 
central component (see Forman et al. 2001a) demonstrates that the central X-ray peak is 
extended and not due to the presence of a nuclear source, thus solving any ambiguity present 
in the HRI data. 

Figure 15 confirms the presence of a clump of hot gas in the core of the companion 
galaxy NGC 508. The count rate measured with Chandra (3.2 ± 0.6 x 10 -3 cnts s -1 ) is 
consistent with the HRI detection and yields an X-ray luminosity 8 Lx ~ 1.5 ± 0.5 x 10 40 
erg s _1 . In the L Y -L S diagram this places NGC 508 (optical luminosity L B = 1.7 x 1O 1O L 0 
(de Vaucouleurs et al. 1991)) in the range occupied by compact early-type galaxies, whose 
emission includes a large contribution from X-ray binaries. 

The preliminary spectral analysis performed by Forman et al. (2001b) shows that the 
mean halo temperature is close to 1 keV, in agreement with the previously published ROSAT 
PSPC (Kim & Fabbiano 1995) and ASCA data (Matsumoto et al. 1997). Taking advantage 
of the Chandra angular resolution, we further analysed the spectral properties of the central 
component. We extracted the 0.3-10 keV spectrum from a 10 arcsec circle centered on the 
central peak, using the surrounding 10-20 arcsec annulus to extract the background. The 
resulting counts were rebinned so to have at least 15 counts per bin and fitted with the CIAO 
Sherpa package. We tried both a single temperature Raymond-Smith (Raymond & Smith 
1977, and additions) and a MekaL (Mewe, Gronenschild & van den Oord 1985; Mewe, Lemen 
& van den Oord 1986; Kaastra 1992; Liedahl, Osterheld & Goldstein 1995) spectral model 
obtaining the results reported in Table 6. The mean gas temperature is kT ~ 0.8 keV for 
the Raymond-Smith model and ~ 0.7 keV for the MekaL one (Figure 16), with the MekaL 
model yielding slightly better fits than the Raymond-Smith one. The metal abundances 
are subsolar but compatible with solar within 3o\ The absorbing columns (nH) are 2-3 
times higher than the galactic one (5.3 x 10 20 cm -2 , ) even though with large uncertainties. 
Fixing the abundance to the solar one yields slightly higher temperatures and lower absorbing 
columns, while fixing the nH to galactic values produces metal abundances higher than solar. 
However, even though the small statistics produces large errors, the plasma temperature is 
quite well constrained as shown by the confidence contours in Figure 16 (right panel). 

We further tried to reduce the extraction region to a 5 arcsec circle, to find out if there 
are gradients in the spectral parameters of the central component. The results of these fits 


8 The error includes the uncertainty on the assumed spectral model, ranging from a thermal plasma with 
variable metallicity to a 5 keV Bremsstrahlung spectrum 
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are shown in the right section of Table 6. The plasma temperatures are slightly lower than 
those extracted within 10 arcsec for all models while the value of the other parameters are 
very uncertain due to the small statistics. However the results of the fits suggest that the 
metallicity may be higher in the galaxy center. 

We are aware that to correctly estimate the hot plasma properties we should perform 
a more complex analysis, e.g. a full spectral deprojection of the different halo components, 
since there is a temperature gradient in the NGC 507 halo. Several authors have shown 
that an incorrect modelling of the X-ray emission (e.g. single instead of multi-temperature 
models) leads to systematic errors in the gas parameter estimates (Buote 1999; Matsushita, 
Ohashi & Makishima 2000). However in our case the small number of events does not allow a 
more complex analysis of the nuclear region. We tried to add the contribution of unresolved 
discrete sources in the form of a 5 keV Bremsstrahlung or a power-law component but the 
fit yields a negligible contribution and, however, does not significantly affect the best fit 
temperature. The best-fit temperature does not change even if we restrict the fit to the 
0.3-4 keV range, where the plasma emission is the dominant contribution. In conclusion we 
believe that the central component temperature is well constrained to be ~ 0.7 keV. 


3. DISCUSSION 

3.1. Origin of the Central Component 

The surface brightness distribution shown in Figure 3 revealed the presence of a bright 
X-ray emission peak centered on the optical galaxy. The surface brightness profiles (Figure 
7 and 8) show that the peak dominates the X-ray emission within the central 10 arcsec (3 
kpc). The HRI data alone do not allow to distinguish between a pointlike or extended nature 
of the central peak, since the comparison of the X-ray profile with the HRI PRF (Figure 
7) and the small core radius of the central component obtained from the bidimensional fit 
(Table 2) are marginally consistent with the presence of a nuclear X-ray source. However 
the Chandra surface brightness profile obtained by Forman et al. (2001a) demonstrates the 
extended nature of this component. 

The comparison between the optical and X-ray profile (Figure 8) seems to rule out the 
possibility that the central peak is due to the unresolved emission of galactic X-ray binaries. 
In this case the X-ray profile of the central component would be close to the optical one, 
since the X-ray source population should follow the stellar distribution. The higher Chandra 
angular resolution and the spectral fit of the central component (§ 2.8) further showed that 
the emission is mainly due to hot plasma and that the contribution of unresolved point 
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sources is small, supporting the gaseous nature of the central peak. 

The presence of a central X-ray excess in galaxies and clusters has been usually asso¬ 
ciated to the presence of a cooling-flow. Kim & Fabbiano (1995) noticed that significant 
mass deposition should be present within the cooling radius since the cooling times of the 
gaseous halo are smaller than the Hubble time (Figure 12). Their estimate, rescaled to our 
adopted distance, yields a mass inflow rate of ~ 20 M 0 yr _1 . If the central X-ray peak was 
due to the gas deposited by the cooling flow, we would expect to find ~2x 10 11 M 0 within 
250 arcsec from the galaxy center while the gaseous mass of the central component is only 
~2x 10 9 (Figure 13, left panel) within the same radius. This suggests that the cooling flow 
is inhomogeneous (Kim & Fabbiano 1995) so that a large amount of cooling gas is deposited 
at large radii and does not reach the halo center. However, if we assume an inhomogeneous 
flow where M oc r (Fabian 1994) and extrapolate the mass deposition down to r = 10”, 
where the contribution of the central component dominates, we still obtain a deposited mass 
of ~ 10 10 M q while we observe less than 10 8 M 0 . This result holds also if we take into 
account all the uncertainties included in the calculation such as the gas filling factor and 
the age of the cooling flow. Moreover we must consider that stellar mass losses are expected 
to have injected an additional 10 9 M 0 of gas in the central 10” (see Figure 13, left panel), 
further increasing the disagreement. 

Thus, if the cooling flow model is correct, most of the deposited gas must have cooled 
out of the X-ray emitting phase. The fate of the cooling gas has been largely debated in 
literature (see for instance Fabbiano 1989; Fabian 1994; Rangarajan et ah 1995; Fabian 2001) 
and is not clear yet. In NGC 507 the column density profile derived by Kim & Fabbiano 
(1995) is consistent with the galactic value and shows that this mass is not in the form 
of cold absorbing gas. One possibility is suggested by Buote (2000b), who found that a 
large amount of warm (kT ~ 10 5-6 keV) ionised gas may be present in the halo of NGC 
507, consistent with the one expected to be deposited by the cooling flow. In this case we 
must address the question of what prevents this gas to rapidly cool to lower temperatures. 
Alternatively a significant amount of gas may have been removed by star formation. This 
hypothesis however, is challenged by the recent failure to observe the cooling gas in the 
center of galaxies and clusters supposed to host large cooling flows (Matsushita et ah 2002; 
Peterson et ah 2001; Kaastra et al. 2001; Tamura et ah 2001). 

The problems mentioned above have led many authors to stress the failure of the stan¬ 
dard cooling flow models and the need to re-analyse the nature of the central X-ray excess 
(e.g. Ikebe et al. 1999; Makishima et al. 2001). The multi-component bidimensional fit of 
the X-ray surface brightness distribution performed in the present work for NGC 507 and 
by Paolillo et al. (2002) for NGC 1399 has shown that the central component is centered on 
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the optical galaxy even though the outer halo is usually asymmetric and displaced several 
kpc (7 kpc for NGC 507, 5 kpc for NGC 1399) from the optical centroid. This suggests that 
the distribution of the denser central gas is related to the underlying stellar distribution. 

Following Makishima et al. (2001) and Matsushita et al. (2002) we compared the total 
over stellar mass profiles of several X-ray bright Early-Type galaxies extracted from liter¬ 
ature. Figure 17 shows that the X-ray halos are characterised by two different dynamical 
ranges: a central region, within 2-3 r e , where the total over stellar mass is almost constant 
with 1 < M Tot /M+ < 3, and an outer region where the M Tot /M* ratio increases steeply 
indicating the presence of massive dark halos. The flat central region correspond to the one 
dominated by the X-ray excess. The Mr 0 t / M* ratio close to 1 indicate that the dark matter 
contribution in the galaxy core is usually small; however even when M Tat /M k is higher, the 
constant value for r < 3r e suggests that the total mass is still related to the underlying 
stellar distribution. The correlation of the ‘knee’ in the Mr at / M* curves with the effective 
radius further supports the association of the central X-ray peak with the optical galaxy. In 
the case of NGC 507 this trend is not as evident as for the other galaxies of the sample for 
the lack of hydrostatic equilibrium in the halo core discussed in § 2.6. However we notice 
that for r > r e , where r e is reported in Table 1) the M Tat /M*. ratio is very similar to those 
of the other galaxies. 

A second point which was usually invoked as a proof of the cooling flow scenario was the 

temperature drop in the halo center. However Brighenti & Mathews (1997) argue that the 

temperature drop within 4 effective radii (r e ) observed in many galaxies, including NGC 507 

and NGC 1399, is not a natural result of cooling flow models. In a recent study of a sample 

of early-type galaxies, Matsushita (2001) explained the properties of X-ray compact galaxies 

as due to the gas injected into the ISM by stellar mass losses. The X-ray luminosities and 

temperatures are compatible with the expectations obtained considering stellar motions as 

2 

the main source of heat, leading to a mean value of /3 spec = r ~ 1. X-ray extended 
galaxies, instead, have a mean halo temperature which exceeds the value expected if the 
gaseous component is in equilibrium with the stellar population (/3 spec ~ 0.5). This difference 
between compact and extended galaxies is reduced if we restrict the analysis to the region 
within ~ 1 r e from the galaxy center, due to the temperature gradient observed in most 
X-ray extended galaxies of the sample. NGC 1399 has the f3 spec value closest to unity 
(/ 3 spec ~ 0.9) among those included in the Matsushita (2001) sample. This fact is easily 
explained considering the analysis performed by Paolillo et al. (2002) which shows that the 
central component dominates the X-ray emission within 1.5r e . 

In the case of NGC 507 (which is not in the Matsushita (2001) sample), if the gas is in 
equilibrium with the stellar component, we expect to find a central temperature close to 0.7 
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keV (assuming a = 330 km s _1 , Prugniel & Simien (1996)), which is in excellent agreement 
with the spectral determinations based on the Chandra data (§ 2.8). The lower central 
temperatures would thus be explained by kinematical heating from stellar mass losses rather 
than by the cooling flow model. In this scenario it is easy to explain why X-ray extended 
galaxies have /3 spec < 1 since to correctly sample the galaxy core we need an angular resolution 
high enough to isolate the temperature of the central component from the one of the outer 
halo, which is usually higher. 

To further confirm the association of the central component with the stellar population 
we computed the X-ray luminosity, within 4 r e , from the central component of NGC 507 and 
NGC 1399 and plotted the result (Figure 18) on the L x {< 4 r e ) vs L B a 2 diagram shown in 
Figure 4 of Matsushita (2001). The Figure shows that, removing the contribution of the 
external halo, the luminosity of the central component falls in the range occupied by X-ray 
compact galaxies whose properties, according to Matsushita (2001), are well explained by 
kinematical heating of the gas supplied by stellar mass loss. 

In conclusion the characteristics of the central X-ray peak analysed in this work for 
NGC 507 and in Paolillo et al. (2002) for NGC 1399 are compatible with the scenario where 
the ISM in the galaxy core is produced by stellar mass losses and is in equilibrium with the 
stellar component. Both gas and stars feel the same potential, mainly due to stellar mass, 
since the amount of dark matter is small. The stellar origin of the gas is further supported 
by the solar metallicity found by Buote (2000a) in the center of several galaxies and groups. 

The scenario that is emerging about the nature of the X-ray excess in galaxies solves 
several problems (origin of the central excess, lower central temperatures, subsolar abun¬ 
dances) largely debated in literature in the last 10 years; however it also reveals our still 
poor understanding of the cooling mechanism: if we assume that the gas remains in the 
hot phase, the X-ray data are inconsistent with the deposition rates predicted not only by 
homogeneous models (a well established fact) but by the inhomogeneous scenario as well. 
Thus either the hot gas must cool out of the X-ray emitting phase, in which case Chandra 
and XMM observations are showing that we are still far from understanding the details of 
the cooling process and the fate of the cooling gas, or there must be some mechanism which 
prevents the hot gas to flow inward. The latter explanation may hold for NGC 507 as we 
will discuss in 8 3.3. 
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3.2. Large-scale Halo structure 

The large-scale properties of the X-ray halo are very different from those of the halo 
core discussed in the previous section. We have shown that the HRI data reveal (§ 2.2) the 
presence of X-ray emission out to 500” (160 kpc) from the galaxy center (Figure 10) and more 
extended toward NGC 499. The mass profiles (Figure 13) show that for r > 80” (25 kpc) 
the total mass greatly exceeds the visible one. While for some galaxies (e.g. Paolillo et al. 
2002) it can be argued that the hydrostatic equilibrium assumption (on which the total mass 
estimates are based) is not valid due to interactions with the surrounding environment, the 
good agreement of the NGC 507 ratio with those derived from literature (Figure 

17) shows that the prevalence of dark matter outside a few effective radii is a common feature 
of many early-type galaxies. 

If the hot gas distribution traces the gravitational potential this dark matter is dis¬ 
tributed preferentially on group scales rather than being associated with the dominant galaxy 
of the cluster. In fact the study of the X-ray surface brightness distribution showed that 
the optical galaxy is displaced with respect to the gaseous halo. This result is confirmed 
independently by the bidimensional fit, which shows that the extended component centroid 
is located ~ 22” (7 kpc) South-West of the optical galaxy (Figure 6, left panel), and by 
the wavelets algorithm, which centered source No. 10 (associated with the NGC 507 halo) 
South-West of source No.9 (coincident with the central peak, see left panel of Figure 3). 
A very similar situation, where the central X-ray peak is centered on the optical galaxy 
while the external halo is displaced from the optical distribution, was found by Paolillo et 
al. (2002) for NGC 1399, thus supporting the idea that the large-scale distribution traces 
the gravitational potential of the galaxy group. In the case of NGC 507 the displacement 
toward NGC 499 seems compatible with the distribution of galaxies within the NGC 507 
group studied by Barton, de Carvalho & Gcllcr (1998). 

The comparison of the halo temperature outside a few effective radii and the stellar 
velocity dispersion in several early-type galaxies yields values of (3 spec ~ 0.5 (Matsumoto 
et al. 1997; Brown & Bregman 1998, 2000; Matsushita 2001) thus supporting the need for 
additional heat to the one provided by stellar mass losses alone. Sources of heat such 
as a central active nucleus or supernovae explosions, seem not enough to account for the 
observational data (e.g. Matsumoto et al. 1997) and would however affect in first place 
the halo center. The alternative explanation proposed by Makishima et al. (2001) and 
Matsushita (2001) is that X-ray extended galaxies are at the bottom of large potential 
structures corresponding to galaxy groups, which provide the additional heat by gravitational 
infall. This scenario is consistent with our analysis which shows that NGC 507 and most 
X-ray bright early-type galaxies, lies at the center of extended dark halos. 
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The displacement of NGC 507 with respect to the large-scale gas distribution may be 
produced by the motion of the stellar body within the extended gaseous halo. The central 
component would then be centered on the optical galaxy since it is due to gas of stellar 
origin in equilibrium with the stellar component (§ 3.1). The tails seen in the adaptively 
smoothed images (Figure 4 and Figure 15) may represent cooling wakes produced by such 
motion, similar those found by David et al. (1994) or Merrificld (1998). Forman et al. (2001b) 
suggested a similar origin for the sharp Southern and Eastern edges in the gas distribution. 
The smaller North-Western tail departing from the central component found in the surface 
brightness map (Figure 3) may be produced by such motion as well, even though we can not 
exclude that it simply represents a layer of gas compressed by the radio plasma on the edge 
of the Western radio lobe. 

Alternatively tidal interactions between NGC 507 and NGC 499 may be responsible for 
the asymmetries in the hot gas distribution. Such interactions may produce large density 
fluctuations and the formation of filamentary structures (D’Ercole, Recchi & Ciotti 2000). 
The presence of ‘clumps’ of cooler gas in the halo, suggested by Kim & Fabbiano (1995) (see 
§ 3.4) may support such scenario. 

These considerations do not necessarily exclude each other since a combination of dif¬ 
ferent effects is likely to be present in high density environments. 


3.3. X-ray/Radio Correlations 

We showed in § 2.2 that the central NGC 507 halo has a very complex morphology. 
Figure 3 revealed the presence of secondary emission peaks and cavities in the X-ray emis¬ 
sion. The weak dependence of the HRI spectral response on the temperature near 1 keV 
(Clarke, Harris & Carilli 1997), suggests a correlation between surface brightness and density 
fluctuations of the hot gas. The comparison of these X-ray features with the radio maps of 
Figure 5 supports the scenario where the hot gas has been displaced by the pressure of the 
radio-emitting plasma. In fact the very good agreement between the shape of the Western 
radio lobe and the X-ray cavity indicates that the lobe may have produced a low-density 
zone in the gaseous medium, compressing the gas at the edge of the lobe (see Clarke, Harris 
& Carilli 1997). This interaction is also likely to be responsible for the N-S elongation of the 
lobe since the radio plasma expanded avoiding the denser Western and Southern regions. 

The presence of such interactions is not so evident in the Eastern lobe. This may be 
due to projection effects and/or to the fact that the Eastern lobe seems to expand in a less 
dense environment (the X-ray centroid of the extended component is displaced by ~ 20” =6 
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kpc South-West of the nuclear radio source). The Radio/X-ray comparison in left panel 
of Figure 5 shows that a secondary emission peak lies in the center of the Eastern radio 
lobe. We followed Feigelson et al. (1995) to calculate the expected X-ray emission due to 
Inverse Compton scattering of Cosmic Background photons into the X-ray band, comparing 
the minimum energy held strength Bme with the magnetic held Bjc required to produce a 
measured ratio of radio to X-ray hux. We used the radio hux for the Eastern lobe measured 
by de Ruiter et al. (1986) and we measured the X-ray hux density extracting counts from 
a 25” circle centered on the secondary emission peak. We found that the expected X-ray 
emission falls an order of magnitude below the measured excess for any choice of the involved 
parameters, and thus Inverse Compton scattering is unlikely to produce the emission peak. 
Conversely, the higher resolution radio contours shown in right panel of Figure 5 suggest 
that the bulk of the radio emission may fall south of the peak, in a low emitting region. 
In this case both the secondary emission peak and the SE ‘tongue’ found in the adaptively 
smoothed image (Figure 4) may be produced by hot gas compressed by the radio-emitting 
plasma, similarly to what seen for the Western lobe. 

The thermal confinement of the radio lobes is supported by the interactions between the 
radio and X-ray emitting plasma discussed above. Morganti et al. (1988) studied the pressure 
balance between the radio lobes and the surrounding ISM, finding that the thermal pressure 
greatly exceeded the radio one (Pthermai / Pradio >7). However their estimate was based on 
Einstein data, which resulted in a temperature twice the ROSAT PSPC value. Taking into 
account the ROSAT temperature reduces the disagreement by more than 50%; the remaining 
excess may be explained by the departures from the minimum pressure conditions discussed 
in their work. 

If the secondary peaks are due to hot gas compressed by the radio emitting plasma we 
can estimate the mass and energy contained in such overdensities. Assuming that the excess 
emission is clue to spherical bubbles with radii of 25” (8 kpc) and 17” (5 kpc) for the western 
and eastern lobe respectively we calculate a mean density of ~ 6 — 7 x 10 “ 3 cm -3 and a 
total mass of a few 10 8 M 0 . The lack of energy resolution of the HRI data does not allow us 
to estimate the temperatures of the overdensities; however the Chandra temperature map 
obtained by (Forman et al. 2001b) suggests that the bubble temperature is not significantly 
different from the one of the surrounding medium. Assuming kT ~ 1 keV we obtain a mean 
pressure P ~ 2 x ICE 11 erg cm -3 . Thus the density and pressure of the ISM contained 
in the bubbles are approximately twice those of the surrounding gas. An estimate of the 
energy required to produce such peaks can be obtained from the excess of thermal energy 
3 nkTV contained in the bubbles with respect to the ambient medium. We obtain ~ 10 5 ' ergs 
summing both bubbles. This estimate is consistent within a factor of ~ 2 with the energy 
required to produce the western cavity 7/(7 — 1 )PV, which is a good agreement considering 
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the uncertainties involved in the calculation. Since the bubbles are overpressurised they will 
expand in the surrounding ISM, releasing part of their internal energy as mechanical work. 
Assuming that the expansion is adiabatic all the excess internal energy will be transferred 
to the ISM. The energy injected by the central radio source offers a possible explanation to 
the inconsistencies of the cooling flow scenario discussed in § 3.1: in fact if radio sources are 
intermittent phenomena with evolutionary timescales of 10' yr we obtain an energy input 
in the ISM of ~ 10 42 erg/s, which is comparable with the energy radiated by the X-ray 
emitting gas within the central 1.5 arcmin of the NGC 507 halo. This energy input should 
significantly affect the gas cooling time preventing the hot gas to inflow and deposit in the 
halo core. This result is in excellent agreement with the simulations performed by Rizza et 
al. (2000). In their models the interaction of the radio lobe with the ISM produces X-ray 
excess emission around the radio lobes of ~ 20 — 30% such as those observed in NGC 507. 
Moreover they estimate that the kinetic energy of the radio jets is comparable or higher than 
the energy lost by the cooling flow via thermal emission. 

Finally we notice that Canosa et al. (1999), using the ROSAT HRI data, found an upper 
limit on the X-ray flux of a possible nuclear source that exceeds, by more than two orders of 
magnitude, the emission expected from the Radio data. The extended nature of the central 
peak implies that the flux of the nuclear source is much smaller than thought before (Forman 
et al. 2001a), thus reducing the disagreement between the radio and X-ray data. 


3.4. Discrete Sources 

Kim & Fabbiano (1995) revealed the presence of many additional pointlikc sources in 
the NGC 507 halo. The co-added spectrum of these sources was found to be well fitted by a 
soft thermal model (kT=0.52) with low abundance (0.0-0.2), leading the authors to speculate 
that they could be cooling clumps in the gaseous halo. As discussed in § 2.7 we used the 
wavelets algorithm to investigate the nature of discrete sources in the HRI held (Figure 14). 
After excluding NGC 507 (source No.10), NGC 499 (No.11), we cross-correlated the HRI 
sources (Table 3) with those detected in the PSPC data (Table 2 of Kim & Fabbiano (1995)). 
Only 5 out of the 17 PSPC sources sources falling within the HRI held 9 are detected in the 
HRI data. We compared the fluxes measured by the two instruments, since any temporal 
variability would exclude the extended nature of these sources. As shown in Tables 4 and 5 
two sources varied signihcantly between 1993 and 1995. One more object (Source No.7) is 
likely to be associated with a background radio source. 


9 Sources No.1,2,19,21 of Kim & Fabbiano (1995) fall outside or on the edge of the HRI FOV. 
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In conclusion ~ 20% of the sources detected by Kim & Fabbiano (1995) are likely to 
be compact objects; this number must be considered as a lower limit since 4 additional 
HRI sources possess an optical counterpart candidate. This result is not necessarily in 
disagreement with the soft spectrum reported by Kim & Fabbiano (1995) since Sarazin, 
Irwin & Bregman (2001) claimed the discovery of supersoft sources in NGC 4697 (even 
though their sources are preferentially located near the galaxy center). The nature of the 
remaining sources is still uncertain. If the emission is due to accreting binaries associated 
with NGC 507 these systems must host massive black holes since their luminosity would 
exceed, by at least an order of magnitude, the Eddinghton luminosity for a 1.4 M 0 Neutron 
Star. In such case they may be placed in globular clusters, given the large distances from the 
galaxy center (see for instance Angelini, Loewenstein & Mushotzky 2001). However, we can 
not exclude the diffuse origin of many sources since the irregular structure of the gaseous halo 
and the filamentary features seen in the halo core suggest that significant density fluctuations 
may be present in the hot gas distribution, similar those found for NGC 1399 (Paolillo et al. 
2002 ). 


4. Conclusions 

The ROSAT HRI data showed, in agreement with previous studies, a bright X-ray halo 
surrounding NGC 507. The large-scale surface brightness distribution is irregular and more 
extended in the North direction. The halo core revealed a complex morphology: there is a 
main X-ray peak, coincident with the position of the optical galaxy, and several secondary 
peaks, due to the interaction between the radio lobes and the surrounding ISM. 

Modelling the halo surface brightness with a bidimensional double Beta model, we found 
that the extended halo centroid is displaced ~ 22” (7 kpc) South-West of the optical galaxy. 
This result suggests a different origin for the central and the extended components. The 
ROSAT and Chandra data show that the central X-ray peak is produced by dense hot 
gas in the galaxy core. Even though the cooling time of the halo within 250” (80 kpc) is 
shorter than the Hubble time, the mass of the gas present in the central component falls 
two orders of magnitude below the amount expected to be deposited by both homogeneous 
or inhomogeneous cooling-flow models. The core luminosity and temperature are instead in 
good agreement with the scenario where the central X-ray peak is due to gas of stellar origin 
kinetically heated by stellar mass losses. This result is supported by the solar metallicity 
found by Buote (2000a) and suggests that the core of bright galaxies is similar to the halo 
of X-ray compact early-types studied by Matsushita (2001). 

The displacement of the cluster halo from the optical galaxy and the filamentary struc- 



23 


tures observed in the halo core suggest that the galaxy may be slowly moving within the 
group potential. The total mass profile derived from the bidimensional model shows that 
a significant amount of dark matter is required at large radii, if the gas is in hydrostatic 
equilibrium. The dark halo extends on group scales and is likely associated with the whole 
cluster rather than with the optical galaxy. 

Comparing with data in the literature, we find that the separation of the gaseous halo 
in two different dynamical ranges (core associated to the stellar distribution; external halo 
tracing the larger potential structures corresponding to galaxy groups and clusters) is a 
common feature of many X-ray extended galaxies. While this result may solve the problem 
of the different X-ray properties among X-ray extended and compact galaxies, since the 
latters are missing the large-scale component (Ikebe et al. 1999; Makishima et al. 2001; 
Matsushita 2001), it makes more difficult to explain the missing deposited mass expected 
from the short cooling times. 

In NGC 507 we found that, if radio sources are intermittent phenomena with timescales 
of 10' yrs, the amount of energy injected in the ISM is large enough to prevent the gas from 
cooling in the halo core and may explain the failure of the cooling flow scenario. 

Finally we found that ~ 20% of the sources detected by Kim & Fabbiano (1995) in 
the NGC 507 halo are due to discrete sources. The nature of the remaining population is 
unclear. If associated with NGC 507, they can be either accreting binaries hosting a massive 
black hole or density clumps of the X-ray halo. 
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Table 1: The ROSAT HRI observation of the NGC 507/NGC499 field. 


name 

Field center 

NGC 507 

(R.A., Dec) 

01 h 23 m 39 s 33°15’36” 

sequence id. 

RH600680n00 

Exp. time 

28310 s 

obs. date 

1995 Jan 14 

P.I. 

D.-W. Kim 

No. OBI 

2 

Distance 

65.8 Mpc 

r e (b) 

53 arcsec 

N h {c) 

5.3 x 10 20 cm -2 


(“^Estimated assuming H$ = 75 km s” 1 Mpc -1 and a velocity of 4934 km s -1 (Huchra, Vogeley & Geller 1999) 
^effective radius (de Vaucouleurs et al. 1991) 
b)Galactic line-of-sight column density (Stark et al. 1992) 


Table 2: Best-fit parameters for the bidimensional halo model. 


Component 

Center Position 

ro 

P 

e 

9 

xl 

V 


R.A. 

Dec 

(arcsec) 



(rad) 


(d.o.f.) 

Central 

01 h 23 m 39.8 s 

33°15’26” 

0.22±0.02 

0.58±0.01 

O.O^b 

0.0 ( b 

0.5 

73 

Extended 

01 h 23 m 38.7 s 

33°15’07” 

49± 9 

0.48±0.04 

0 . 0 (a) 

0 . 0 (a) 

0.9 

776 


The parameters e and 9 have no error because they were held fixed during the fit (see discussion in text). 
Note - uncertainties are ler confidence level for 5 interesting parameters 






Table 3: Discrete Sources - Results of the Wavelets algorithm 


Source 

No. 

R.A. Dec. 

(J2000) 

3cr Radius 
(arcsec) 

Counts 

Max S/N 

Count Rate 
(10 -4 Cnts sec -1 ) 

rpow.law 

J HRI 

fRS (a) 

Jhri 

1 

1:23:38.6 

+33:21:48 

8.7 

13+6 

4.3 

5+2 

35+14 

20+8 

2 

1:23:07.5 

+33:15:23 

12.2 

31+9 

7.0 

11+3 

76+21 

45+12 

3 

1:23:11.5 

+33:10:58 

12.2 

37+10 

8.0 

13+4 

90+28 

53+16 

4 

1:23:25.0 

+33:25:31 

17.3 

47+13 

7.6 

17+4 

118+28 

69+16 

5 

1:23:43.2 

+33:15:21 

17.3 

33+11 

4.8 

12+4 

83+28 

49+16 

6 

1:22:58.7 

+33:22:43 

24.5 

39+12 

5.2 

14+4 

97+28 

57+16 

7 

1:23:44.6 

+33:26:08 

24.5 

33+11 

4.6 

13+4 

90+28 

49+16 

8 

1:23:43.4 

+33:32:23 

34.6 

46+14 

4.9 

21+6 

145+42 

65+20 

9 

1:23:35.8 

+33:15:03 

34.6 

104+28 

6.1 

37+10 

256+69 

151+41 


(“) Flux in the 0.1-2.4 KeV band in units of 10 15 erg s 1 cm 2 . 


Table 4: Discrete Sources - Comparison with PSPC results and catalogues 


Source 

No. 

PSPC (a) 

No. 

fRS (a,b) 
JPSPC 

fRS (6) 
Jhri 

Time Var. 

long-term^) 

Opt.id. 

NED obj. 

1 

- 

- 

20+8 

- 

n 

- 

2 

8 

21 ± 10 

45+12 

n 

y 

- 

3 

10 

171 + 20 

53+16 

y 

y 

- 

4 

12 

91 ± 12 

69+16 

n 

y 

- 

5 

- 

- 

49+16 

- 

n 

- 

6 

7 

81 ± 12 

57+16 

n 

y 

- 

7 

13 

43+10 

49+16 

n 

y 

NVSS J012345+332554 

8 

- 

- 

65+20 

- 

n 

- 

9 

- 

- 

151+41 

- 

n 

- 


(a) reported in Table 2 of Kim & Fabbiano (1995) 

Flux in the 0.1-2.4 KeV band in units of 10 -15 erg s -1 cm -2 . 

O Difference between PSPC and HRI fluxes: ‘y’ means variable source (> 3cr), ‘n’ means no variability 
detected. 






Table 5: PSPC-HRI Sources Comparison - Upper limits from HRI data 


PSPC(“( 

No. 

fRS (a 
JPSPC 

,6) fRS 
JHRI 

fRS 

JHRI 

3a upper limit 

Time Var. 
long-term( c ( 

3 

79 ± 15 

14 ±25 

90 

- 

4 

36 ± 10 

21 ±24 

93 

- 

5 

61 ± 14 

40 ±24 

111 

- 

6 

300 ± 27 

116 ±42 

243 

y 

9 

45 ±11 

8 ±23 

77 

- 

11 

42 ±11 

14 ±24 

86 

- 

14 

38 ± 10 

32 ±25 

108 

- 

15 

32 ± 10 

9 ±21 

72 

- 

16 

28 ± 10 

9 ±22 

75 

- 

17 

51 ±11 

11 ±22 

78 

- 

18 

46 ± 13 

67 ±23 

135 

- 

20 

60 ± 14 

19 ±21 

83 

- 


(“(As reported in Table 2 of Kim & Fabbiano (1995) 

Flux in the 0.1-2.4 KeV band in units of 1CK 15 erg s _1 cm -2 . 

(“(Difference between PSPC flux and HRI upper limit: ‘y’ means variable source (> 3a), means 
no variability detected. 


Table 6: Best-fit spectral parameters for the central component from Chandra ACIS-S data. 


Model 


(r < 10 

arcsec) 



(r < 5 

arcsec) 



kT 

Abundance 

nH 

Xledl V 

kT 

Abundance 

nH 

xleJv 


(keV) 

(Z/Z Q ) 

(10 20 cm -2 ) 


(keV) 

{Z/Zq) 

(10 20 cm -2 ) 


Raymond-Smith 

0.82 ± 0.03 

0.48™ 

16 ±5 

0.89/60 

0.74 ±0.05 

0.48-;- 

ie±i 

0.74/39 

Raymond-Smith 

0.84 ±0.02 

l(“) 

11 ±4 

0.91/61 

0.77 ±0.04 

l(«) 


0.75/40 

Raymond-Smith 

0.86 ±0.01 

1 q+0.8 

i -°-0.4 

5.3(“( 

0.96/61 

0.81 ±2 

i-9;H 

5.3(°) 

0.81/40 

MekaL 

0.73 ± 0.04 

Q QQ+ 0 - 21 
u - oy —0.12 

21 ±5 

0.78/60 

0.68 ±0.04 

o.stJi 

21 ±6 

0.62/39 

MekaL 

0.77 ±0.04 

i(«) 

17 ±6 

0.83/61 

0.68 ±0.04 

l(«) 

21 ±1 

0.60/40 

MekaL 

0.81 ±0.02 

1 co+1-2 
±.00-0.5 

5.3(°( 

0.94/61 

0.74 ±0.04 

q+ 123 
°-l .4 

5.3(“( 

0.74/40 


(“(Fixed 

Note - Errors are la for all interesting parameters. 
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Fig. 1.— Raw HRI image of the NGC 507/NGC 499 field. The data are rebinned in 5 x 5 
arcsec pixels. North is up and East is to the left. 
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Right Ascension (J2000) 

Fig. 2.— Exposure corrected image of the center of the HRI held. The data are rebinned in 
5x5 arcsec pixels and smoothed with a gaussian filter of a — 15”. The X-ray contours are 
spaced by a factor of 1.2 with the lowest one at 6.4 x 10 -3 cnts arcmin -2 s _1 . Crosses show 
the position of sources detected by the wavelets algorithm (Table 3). The white dashed box 
is the region enlarged in Figures 3, 5 and 6. In the figure North is up and East is to the left. 
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Fig. 3.— Left: 2.5 x 2.5 arcsec/pixcl image of the nuclear region of NGC 507. The image 
is convolved with a gaussian filter of a = 5”. Contours are spaced by a factor 1.3 with the 
lowest one at 1.4 x ICC 2 cnts arcrnin -2 s^ 1 . The dashed circles show the region of maximum 
S/N ratio for the sources detected by the wavelets algorithm. A higher quality color version 
of this figure is available in the electronic edition. Right: X-ray brightness contours overlaid 
on the 1 arcsec/pixcl Digital Sky Survey (DSS) image (logarithmic grayscale). 
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Fig. 4.— Adaptively smoothed image of the center of the NGC 507 halo. The data are 
rebinned in 1 x 1 arcsec pixels. Contours are spaced by a factor 1.2 with the lowest one at 
8.3 x 10~ 3 cnts arcmin -2 s _1 . The dashed box is the region enlarged in Figures 3, 5 and 6. 
A higher quality color version of this figure is available in the electronic edition. 
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Fig. 5.— Left: 13” FWHM radio contours superimposed on the X-ray nuclear region of 
NGC 507. The X-ray image is the same shown in the left panel of Figure 3. Right: Higher 
resolution (3.5” FWHM) radio contours superimposed on the X-ray nuclear region of NGC 
507. Higher quality color versions of these figuree are available in the electronic edition. 
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Right Ascension (J2000) Right Ascension (J2000) 


Fig. 6.— Left: X-ray contours of the bidimensional model superimposed on the NGC 
507 central halo. The model contours are spaced by a factor 1.3 with the lowest one at 
9.6 x 10 -3 cnts arcmin -2 s -1 . Right: X-ray contours superimposed on the residuals of the 
bidimensional model. The contour levels are the same of the right panel of Figure 3. The 
residuals range from —1 x 10 -3 (white) to +2 x 10 -3 cnts arcmin -2 s -1 (black). The crosses 
show the grayscale level corresponding to zero. Higher quality color versions of these hguree 
are available in the electronic edition. 
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Fig. 7.— HRI surface brightness profile of NGC 507 (thin continuous line). The dashed (dot- 
dashed) line represent the SMB background map profile after (before) rescaling to match 
PSPC counts. Filled circles show the background subtracted profile obtained using the 
rescaled SMB background. The PRF profile and the uncertainties due to residual aspect 
errors, are shown respectively as the thick continuous line and the shaded region. 
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Fig. 8.— Radial profile of the bidimensional HRI model (continuous line). The central and 
extended components are represented respectively by the dashed and dot-dashed lines. HRI 
(PSPC) counts are shown as filled circles (empty squares). The dotted line shows the V 
band surface brightness profile measured by Gonzalez-Serrano & Carballo (2000). 
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Fig. 9.— Comparison of the azimuthally averaged (filled circles, continuous line) X-ray 
surface brightness profile (same as the one shown in Figure 8) with the one extracted in the 
NE quadrant (empty circles, dashed line). 
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r (arcsec) 


Fig. 10.— NE (20° <P.A.< 110°) and SW (110° <P.A.< 200°) radial surface brightness 
profile of the NGC 507 X-ray halo. 







41 



Fig. 11.— Deprojected density profile of the NGC 507 gaseous halo (continuous line). The 
dashed and dot-dashed lines show, respectively, the contribution of the central and extended 
component. Crosses represent the density obtained using the NE quadrant profile shown in 
Figure 9. 
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Fig. 12.— Cooling time profiles of NGC 507. Symbols have the same meaning as in Figure 
11. The cooling times derived by Kim & Fabbiano (1995) are shown as filled squares. 
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Fig. 13.— Left: Integrated gaseous mass within r. Symbols are the same as in Figure 11. The 
dotted line represents the mass injected in the ISM by stellar mass losses. Right: Integrated 
total mass within r. The continuous lines represent the gravitating mass estimates obtained 
assuming the different temperature profiles described in the text. The total mass estimate 
of Kim & Fabbiano (1995) is shown as filled squares. The dashed line shows the total mass 
obtained using the surface brightness profile extracted in NE quadrant (see Figure 9). The 
shaded (dotted) region shows the stellar mass contribution estimated (extrapolated) from 
the optical profile measured by Gonzalez-Serrano & Carballo (2000), assuming a M/L ratio 
ranging from 6 to 8 M Q /L Q . 
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Fig. 14.— Discrete sources detected by the wavelets algorithm superimposed on the adap¬ 
tively smoothed HRI field. The circles show the region of maximum S/N ratio for each 
source. The position of sources No.5 and 9 is shown in greater detail in the left panel of 
Figure 3. 
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Fig. 15.— Chandra ACIS-S adaptively smoothed image of NGC 507. The black box shows 
the edges of the ACIS-S3 chip. Contours are spaced by a factor 1.2 with the lower one at 
2.4 x 10~ 2 photons arcmin -2 s -1 . A higher quality color version of this figure is available in 
the electronic edition. 
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Fig. 16.— Left: Best-fit MekaL spectral model of the central 10 arcsec region. Residuals 
are shown in the bottom panel. Right: 1, 2 and 3<r confidence levels for the best-fit MekaL 
model. 
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Fig. 17.— Total over stellar mass profiles derived from literature. Total mass estimates are 
derived from X-ray measurements while the stellar contribution is obtained from optical and 
dynamical estimates. 
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Fig. 18.— Lx(< 4 r e ) vs for a sample of Early-Type galaxies. The Figure is identical 

to Figure 4 of Matsushita (2001) except for the addition of NGC 507. The stars represent 
the drop in X-ray luminosity for NGC 507 and NGC 1399 considering only the contribution 
of the central component. The continuous line is the best fit to X-ray compact galaxies while 
the dotted line represents the kinetic heating rate by stellar mass loss (see Matsushita 2001 
for details) 




